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Abstract 
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The results of analytical evaluation of 0{a s ) QCD contributions due to the nonvanish- 
ing quark masses to Tu_^ qf - q - f are presented. The "triangle anomaly" type contributions 
are included. As a byproduct the 0(a^) logarithmic contributions are evaluated. The 
\ results are presented both in terms of running and pole quark masses. The partial de- 

cay modes H — > 66 and H — > cc are considered. The calculated corrections decrease 
Qh. the absolute value of large and negative 0{o? s ) massless limit coefficient by < 1% in 

the intermediate mass region and by 1%— 20% in the low mass region which, however, 
is experimentally ruled out. The results are relevant for H — > tt decay mode for the 
higher Higgs mass region where the mass effects are large and important. The high order 
corrections remove a very large discrepancy between the results for T}j^q f q f in terms of 
running and pole quark masses almost completely and reduce the scale dependence from 
about 40% to nearly 5%. The remaining theoretical uncertainties are discussed. 



1. Introduction The discovery of the Higgs boson, which would be a crucial confirmation 
of the Standard Model (SM), is one of the primary goals of modern high energy physics 
The unsuccessful search of the Higgs particle (H) at LEP 100 has established the present 
lower bound on its mass: M H > 63.5 GeV (95% c.L). On the other hand, the theory is able 
to provide only a broad range of the allowed Higgs mass with the upper limit about 9M\y — 
(87TV2/3GF) 1 / 2 ~ 1 TeV (see, e.g., 0). However, if the standard spontaneous symmetry 
breaking mechanism is correct, then the SM Higgs particle is expected to show up at LEP 200 
or/and LHC (for the review of Higgs search at LEP see, e.g., 0) and possibly at the upgraded 
Fermilab Tevatron 0. 

Prior to the experimental discovery, important constraints on the Higgs properties and on 
its mass can be obtained when the loop radiative corrections are taken into account. Moreover, 
in the case of discovery of the Higgs boson, the precise evaluation of the characteristics of its 
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production and decay processes from the first principles of the theory would be important in 
order to distinguish the SM Higgs from its extended theoretical versions. (For the analysis of 
the loop radiative effects in the SM Higgs phenomenology see 0.) 

The explicit calculations show that the high order radiative effects are large and important. 
For example, the calculated QCD corrections at Mu ~ 120 GeV reduce the Born result for 
H — ► cc by about 20% and for the bb mode by about 42% (see, e.g., and references therein). 
There are similar size effects for H — > gg and even larger for pp — > H + X. 

The leading order QCD correction to the decay width H — ► qfqf, where q$ is the quark with 
flavor / and mass m/, has been calculated by several groups In || this correction has 

been evaluated in the zero quark mass limit. The 0{a^) QCD correction to the T H ^ qf g f has 
been computed in || (see also |H|) in the zero quark mass limit and the results were presented 



in terms of running quark mass. 

As was mentioned, although the QCD corrections depend on Mr, their impact on the 
decay width is very large for the wide range of the Higgs mass. Indeed, the massless limit 
(rrif/MH — > 0) three-loop correction [I(J is about 10% of the Born result and 30% of the 
leading order QCD correction (with a s (Mn) ~ 0.15). On the other hand, presently there 
are no analytical methods to carry out the complete three-loop calculations involving massive 



particles (as an exceptions see JIT], [12]). However, the contributions due to the nonvanishing 
quark masses might be important for modeling the behavior near the flavor thresholds including 
the it threshold. Moreover, the possibility exists that the large 0(a 2 ) QCD correction || fO 



is an artifact of the limit mf — > 0. 

In the present work the 0(a 2 s ) corrections ~ m 2 /Mfj (/ = u, d, s, c, b, t) to the decay width 
rH-, qf q f are calculated. The results are obtained both in terms of running and pole quark 
masses. The calculated ~ rnj/M^ corrections are not entirely negligible in the intermediate 
and low Higgs mass regions. For instance, at M# = (100, 40, 20)GeV the 0(a 2 s ) mass correc- 
tions increase the previous result for T H ^ b i by about (1%,5%,23%). However, the low Higgs 
mass region is experimentally ruled out |J. On the other hand, the mass corrections for the 
it decay mode in the high Higgs mass region are very large. 

2. Preliminary Relations. The standard SU(2) x U(l) Lagrangian density of fermion- 
Higgs interaction is: 

L = -9 Y q f q f H = -{V2G F fl 2 m f q f q f H = -{^/2G F ) l ' 2 jfH (1) 

The two-point correlation function of the scalar currents jf = mfq^qj has the following form: 

n(Q 2 = -s,m f ) = i J e iqx < Tj f (x)j f (0) > d A x (2) 
The decay width can be expressed as the imaginary part of U(s + iO, mf) in the standard way: 



(3) 
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The total hadronic decay width will be the sum over all participating (depending on Mh) 
quark flavors: 

f=u,d,s,... 

For the experimentally accessible quantity - the branching fraction of the particular decay 
mode one has: 

BR(H -> q f q f ) = V H ^ qf - q} / iY^ hadrons + V H ^) (5) 

The full 0(a s ) analytical result for the decay rate of H — > qjq^ in terms of pole quark masses 
looks like §-0: 



rv 



3\[2G F M H 



where: 



8tt 
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m}[ 1 - 



4m 2 x - 



i + 



71 



(6) 



a(P) 3 + 3A/3 2 - 13/3 4 , 21/3 2 - 3 

— log7 + 







16/3 3 



l(3 2 



a{(3) = {\ + (3 2 ) 



4Lz 2 ( 7 - 1 ) +2Lz 2 (- 7 - 1 ) - log 7 log- 



:i+/5) 3 

7 = (1 + - 0), = (1 - 4m 2 /il4) 1 / 2 

and the Spence function is defined as usual: 
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The expansion of the r.h.s of eq.(^) in a power series in terms of small m 2 /Mjj has the following 
form: 



3y/2G F M H 



H^q f q f 



87T 



1 - 6 



m 



f 



Ml 



+ ... 



a s {M H ) 



7T 



M 2 



3-21og— ^" 



m, 



m 2 M 2 H 



8-24 log 



Ml 

2 
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+ 0(a 2 



(7) 



where the period covers high order terms ~ (rUf/Ma) , k = 2,3. 



3. 0(a 2 ) QCD contribution to TH^ qf q f - The expansion of the full two-point correlation 
function (|2|) in powers of mj/Q 2 in the "deep" Euclidean region (Q 2 ^> m 2 ) has the following 
form: 



m 



Tl{Q 2 ,m f ,m v ) = U^Q 2 ) + '-^U m ,(Q 2 ) + £ ^n m? (Q 2 ) + ... 



m 2 Q 2 



Q 



v=u,d,s,c,b 



Q 1 



The last term in the above expansion is due to the certain topological types of three-loop 
diagrams containing virtual fermionic loop. In fact, in these loops the virtual top quark can 
also appear. This issue will be discussed later. 
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In order to evaluate the coefficient functions on the r.h.s of eq.(^j), it is sufficient to write 
the diagrammatic representation for the II (Q 2 , mf, m^) up to the desired level of perturbation 
theory and apply the appropriate projector. To 0(a 2 ) one has: 

n 2 (o 2 )- 1 ( d \ 2n ( d V fc f n (QW.™« B ) ] (q) 

L W<"W > (2n)\(2k)\\dmf J Urn* J I m^i-n-fe) j mf=m?=0 ^ 

where n, k — 0, 1; n + k < 1 and superscript "B" denotes the bare quantities. 

The one-, two- and some of the typical three-loop diagrams contributing to the IT are 
pictured on the fig. 1. 



fig.l Feynman diagrams contributing to the Y H _, q{1}f to 0(a 2 ). 
The total number of three-loop graphs is 18. 



The dimensional regularization formalism 



13| and the minimal subtraction prescription 



14J in its modified form - MS |]15[ are used. The corresponding one-, two- and three-loop 
diagrams have been evaluated analytically, using the special computer program HEPLoops 
, written on the analytical programming system FORM [ 17 1 . More details and the graph- 



by-graph results will be given elsewhere. The graph-by-graph results were summed up with 
an appropriate symmetry and gauge group weights. In the obtained expression for the IT in 
terms of bare quantities one renormalizes the coupling and the quark mass factors: 
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nit 



a. 



6C f 



47T 
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af = 

OCs\' 

4vr/ 



c f 



4tt V 3 3 ) e 



(11Ca + 18C f -47W)- 



(10) 



9 {c A+ lc F - 



™tnY- 

3 ) £ 



(11) 

where, e = (4 — D)/2 is the parameter of dimensional regularization, D is the dimension 
of space-time; N is the number of participating quark flavors. The eigenvalues of Casimir 
operators for the adjoint and the fundamental representations of SU C (3) gauge group are C A = 
3, Cp = 4/3 correspondingly and T — 1/2. 

The obtained expressions for IT at each order of a s are the polynomials with respect to 1/e 
and log [ijjg/Q 2 : . However, there are no terms like (l/£ n )(log/i|^/<5 2 ) fc - They will appear only 
at higher orders ~ m 2 m^ / Q A and represent the infrared mass logarithms. The poles can be 
removed by an additive renormalization. The imaginary part of each of the IT is finite after the 
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regularization is removed. Here the usual way P, |10[ of the introduction of Adler D-function 
fl8|l (the Il-function differentiated with respect of Q 2 ) which is finite and renormalization group 
(RG) invariant, can be avoided. Instead, one can operate directly with the correlation function 
IT. Namely, one should analytically continue it from Euclidean to Minkowski space and take 
the imaginary part at s = Mfj ( eq.@ ). One obtains the following analytical result for the 
standard QCD with SU C (3) gauge group: 



3V2G F M. 
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( OL s 
+ - 
V 7T 



-m 



57T 

10801 
144 
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V24 



/ 
19 



3Q 

C(2) - ?C(3) + 



19 



106 ^ 

3 g Ml + 4 



log 



1 

3 l 



2 , . 11 
3«3) + - 



9 g M 2 



H 



+ \ log 2 
6 



2 ^MS 

M 2 H 



^MS 

M 2 H 



m 



Ml 



a-. s 



7T 



6 + — 40 + 24 log 



Ml 



+ 



2383 
4 ' 

V 18 



162C(2) - 166C(3) + 371 log + 81 log 



4C(2)-4C(3) + 101og^ + 21og 



Ml 



Ml 



S.C.O n 



12 



(12) 



where the Riemann function £(2) = ix 2 / 6 arose from the analytical continuation of log 3 i^raj Q 2 
terms and ((3) = 1.202056903. The last term in eq.(|12|) represents the contributions from the 
three-loop diagrams in fig.l containing the virtual quark loop and the "triangle anomaly" type 
corrections from the graphs in fig. 2. Their contribution vanishes in the massless quark limit. 



fig. 2 The "triangle anomaly" type diagrams. 

The leading three-loop term ~ m 2 (the massless approximation) coincides with the one 
obtained in || TU], while the three-loop results ~ m^m 2 m 2 are new. 

4. Renormalization group analysis. An observable quantity, in particular the cal- 
culated decay width is invariant under the RG transformations and obeys the homogeneous 
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RG-equation: 



H- 



( f 



,mf,m v ,a. 



(13) 



where the QCD /3-function and the mass anomalous dimension are introduced in the standard 
way: 



= 



-a 



fc>0 



7T 



li 2 dm 



rrif dfM 

The known /5-function coefficients are ||19|| : 



k>0 



k+1 



k+1 



Po 



1/, 2 \ n 1/ „ 38 \ „ 1 /2857 5033 Ar 

- ( 11 N ) , p 1 = — 102 N), p2 = —{ N + 

4 V 3 / 16 V 3 / 64 V 2 18 



325 



N- 



and for the 7 m (a s ) one has f20 
7o = 1, 7i 



1 /202 20 



16 V 3 9 



1249 



/2216 160 



V 27 



140 



+ — C(3) JV-— N 



(14) 
(15) 

(16) 

(17) 



The solution of the RG eg . Qp~3f) can conveniently be written as follows: 



^H-+q,q f = r mj(/i) £ 
0<j<i 



7T 



log J 



Mjj ^ "" 



*« + £ 



m 



v IV1 H 



where the coefficients a^, and are the same as the ones in eq.(|T^). Applying the differen- 
tial operator /i 2 d/ d/i 2 to the both sides of the eq . fll8|) and taking into account the RG-invariance 
of TH-*q f q t and the eqs. (^4]JT5D , one obtains at the 0(a s ): 



a n = 27 a 00 , &n = 47 6 



at the 0(a 2 ): 



'oo • 



a 21 = 27ia 00 + (/3 + 27 )ai , 
a-22 = (Po + 27o)on/2 = (P + 27o)7 a o, 

621 = 471600 + (Po + 47o)6io, 
622 = (Po + 47o)6n/2 = 2(A) + 4 7o )7o&oo, 



(19) 



(20) 



and at the 0(af): 

a 3 i = 2(p + 7o)a 2 o + (Pi + 271) ai + 27 2 a 00 , 
032 = (Po + 7o)a 2 i + (Pi + 2j 1 )a n /2 = (p + 7o)[27ia 00 + (Po + 27 )ai ] + {Pi + 27i)7 a 00 , 
033 = 2(/3 + 70)^22/3 = 270(A) + lo)(Po + 27o)«oo/3, 
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An 



2(A) + 2 7 o)6 20 + (A + 471)610 + 4 72 &, 



00 • 



(/5 + 27o)&2i + (/5i+4 7l )^ 



'32 



(A) + 2 7o ) [4 7 i6 00 + (Po + 4 7o )6i ] + (A + 4 7 i)2 7o 6 o, 

633 = 2(A> + 2 7o )6 22 /3 = 4 7o (A> + 2 7o )(/?o + 4 7o )W3. 

c 3 i = 2(/3 + 2 7o )c 20 (21) 

Since the diagrams with the virtual fermionic loop first appear at O(a^), c o = C\j = c 2 i = 
c 22 = 0. The relations (|T9|j20|) provide a powerful check of the O(of) calculations, while the 
relations ([21]) allow one to evaluate the log-terms at 0(a 3 ), without explicit calculations of the 
corresponding four- loop diagrams. With those relations, the information available at present, 
namely the QCD /3-function, mass anomalous dimension and the 2-point correlation function 
up to the three-loop level, is fully exploited. In fact, the similar relations can be derived for 
the correlation function II. However, the RG-equation for II is not a homogeneous one and 
the anomalous dimension function up to the corresponding order of a s is needed. 
The solution of the RG-equation flT3|) at jij^ = Ml has the following form: 



H->q s q s 
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3y/2G F M H 
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5.66667-40 



7T 



35.93996 - 1.35865JV - 



m 2 AM H ) 



+ 12 £ ml{M H )/M 2 H 

v=u,d,s,c,b 



m}{M H ) 
Ml 



129.72924 - 6.00093ivj 



(22) 



The running coupling is parametrized as follows: 



a s {Ml) 



AlogL 



7T 



+ lo S 2 L ~ Pi l0 § L + &A) - Pi) + 0(L" 4 ), 



(23) 



where L = log(M| r /A|j^). For the running mass one has: 



m/(/ii)/m/(/i 2 ) = 0(a s (/ii))/0(a s (/x 2 )), 



(24) 



where, 
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+- 
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a s (/x) 



7T 



Y7i 
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/?o 2 / 
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Pill 
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Pi 7(A 

A 2 7o 



A) 

A7o 
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a s (/x) 



7T 



(25) 



In the eqs. 



-(Pq) all appropriate quantities are evaluated for the A" active quark flavors. 



A" can be determined according to the scale of M#. At present one usually considers N = 5. 



5. TH-^q } q f in terms of pole quark mass. For the heavy flavor decay mode of the Higgs, 
it is relevant to parametrize the decay rate in terms of pole quark mass (see, e.g., M). Below 
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the result ( p2|) will be rewritten in terms of pole quark mass, assuming that heavy quark is 
not exactly on-shell. 

First, one derives the following general evolution equation for the running coupling to 0(al): 
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(<<) 



(A*) 
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(TV) 
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7T 



7T 
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(26) 
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2 25 
-—(N-n) 
72 v ; 



where the superscript n (N) indicates that the corresponding quantity is evaluated for n 
(N) numbers of participating quark flavors. Conventionally (see, e.g., PH), n (N) is specified 
to be the number of quark flavors with mass < fi (< M). However, the eq . (|26|) is relevant for 
any n < N and arbitrary yU and M , regardless the conventional specification of the number 
of quark flavors. The logm///z terms are due to the "quark threshold" crossing effects and 
the constant coefficients 1/6 = fi^' 1 ^ — Pq, 19/24 = p[ k ~^ — (3^ represent the contributions 
of the quark loop in the /^-function. The sum runs over N — n quark flavors (e.g., I = b if 
n = 4 and N = 5). Note that m/ is the pole mass of the quark with flavor /. For the on-shell 
definition of the quark masses the eq. (|26|) changes - the constant —25/72 should be substituted 
by +7/72. The above equation is derived based on the eq. (|23"D , the QCD matching conditions 



for a s at "quark thresholds" p2| , |23| and the one-loop relation between on-shell and pole quark 
masses. The eq.(|26|) is consistent with QCD matching relation at mj(m^) [^J (see also plj ): 

af'- V/fa/)) = aW(m/(m/)) + (ai N f\m f (m f ))f(C A /9 - 17CV/96)/tt 2 (27) 
Here and below Nj is the number of quark flavors u,d,...,f. Note that the nonlogarithmic 



constant at O(a^) in eq.(p6"D will not contribute in further analysis. 

Next, using the scaling properties of the MS running mass and the eq. (|26|) , one obtains the 
following matching condition: 
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where the constant terms are: 1/12 = 7o(/3cj fc ^ — /3q^)/2, 5/36 = 7j K 1,1 — / -f\ K> and 2/9 
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T^log 
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(k) 



(2? 



(fc-i) 



/3o )■ In general the 5(m/,m//) is the finite contribution of the single virtual 



heavier quark with mass rap, entering when one increases the number of flavors from tV — 1 
to N (one can also consider the particular case my = rrif). From the two-loop on-shell quark 
mass renormalization one has ||25|| : 

(29) 



5{m f , m/0 = -C(2)/3 - 71/144 + (4/3) A(m f >/m f ) 



where 
A(r) 



log 2 r + £(2) — ( logr 
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)L + {r) - (1 -r)(l -r 3 )L_ 



(30) 
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Ldr) 



1 / r , logx 
dx- 



x±l 



L±(r) can be evaluated for different quark mass ratios r numerically (table 1). 

One can relate the MS quark mass m/(m/) to the pole mass rrif using the 0{a 2 s ) on-shell 
results of 



m?'\m f ) = m f [l - 4aW(m f )/to + (16/9 - K f ){aff\m f )M\ 



(31) 



where 



3817 
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(2 + log2)C(2)-^C(3)-^ 
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C(2) + 
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+ i £ A 
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( m i\ 
KrrifJ 



(32) 



The first four terms in If^ represent the QCD contribution with iVy massless quarks, while the 
sum is the correction due to the Nf nonvanishing quark masses. 

Combining the eqs. (pl]j25| ) and the eqs. (p6|)-(|3T|), one obtains the relation between the 
MS quark mass m/(M#) renormalized at Mh and evaluated for the iV-flavor theory and the 
pole quark mass mf. 



mf ] {M H ) =m/{l 



aW(M H ) 



7T 
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- + 70 log — 
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16 
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(iV) 



log— + — (/% ; -7o)log — T 



(33) 



Note that N is specified according to the size of M# and has no correlation with the quark 
mass rnj. Thus, for instance, one can apply the eq.(|33|) to the charm mass mf\MH) evaluated 
for five-flavor theory. 

Substituting eqs. ( PB|J5^jT7|JT^ ) into the eq. , one obtains the general form for the decay 



rate T H ^ qf g f in terms of the pole quark masses: 
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(34) 



The above result confirms the asymptotic form (^) of the two- loop exact result @, while the 
0(o!g) expression is new. 
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26 



m. 



1.46GeV 



The numerical values of A(m//m/) defined in the eq.(^) are given in the table 1. Quark 
masses can be estimated from QCD sum rules: m& = 4.72GeV 

m s /13 

for the 



m 



0.27GeV PS[, m u + m d 



0.02GeV For the upper bound of the top 

is used. On the other hand, 



3C 



mass the latest CDF result m t < 174 +10 + 13 GeV 
lower bound more conservative value m t > 131 GeV given by DO [31| is used. The values of 
A(mi/m t ) are given in the table 1 for both the lower and the upper bounds of the top mass. 



Table 1. The numerical values of A(m//m^) 



f -flavor 


/-flavor 


A (777 1 / TH 


t 


word 


< 0.0001 


t 


s 


0.0025-0.0017 


t 


c 


0.0137-0.0091 


t 


b 


0.0435-0.0291 


t 


t 


0.8587 


b 


u or d 


0.0027 


b 


s 


0.0673 


b 


c 


0.3290 


b 


b 


0.8587 


c 


u or d 


0.0084 


c 


s 


0.2045 


c 


c 


0.8587 


c 


b 


1.9015 


s 


u or d 


0.0454 


s 


s 


0.8587 


s 


c 


2.5767 


s 


b 


4.5605 



In order to achieve better numerical precision within the wider range of M# it is more 
appropriate to use full one- and two-loop results (eq.(|)) and the mass corrected three- loop 
result given by the eq.fl34}). For the dominant decay mode H — ► bb with five participating 
quark flavors one obtains: 



3V2G F M H 
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8tt 
(5 ) (M, 



1 - 



Ami 



+ 



a 



{5) (M H ) 



71 



7X 



-2.23039 + 266.13393 



18.13889- 177.66667 



- 2.66667- 32 



Ml 
M\ 



rag 
Ml 



log 



Ml 



mi 



E 

mi<m b 



{mi 



1 - 



+ 0.08333 - 25 



\m b 



+ 12 J2 



Ml 

9 

m. 



Ami 
Ml 



log 



mv <m b Ml 



Ml 
ml 



(35) 



where is defined in eq.( 



10 



At the present LEP lower bound on Higgs mass Mg ~ 60GeV one has: 



r^= 3 ^ Mg m ;((l-0.037) + a ' 5 ' (Mg) (-7.17 + 0.71) + f^ ) ' M ^V(-95.67 + 3.57)} 

87T L 7T V 7T / J 

(36) 

where the second number in each bracket corresponds to the leading order quark mass correc- 
tion, while the first one is the massless contribution. One can see that the mass corrections 
always decrease the large massless contribution. For the lower values of Mh the relative mass 
corrections increase rapidly. For example, at Mh = 40GeV the O(a^) mass correction decrease 
the large and negative massless coefficient by more then 10% and correspondingly increase the 
decay rate by about 5%. 



Finally, one also needs to include weak and electromagnetic contributions |32], [|. However, 
they will not affect further analysis here. 

The dependence of the partial decay rate T H ^ b i on the mass of the Higgs particle Mh is 
plotted in fig. 3. 



fig. 3 The partial decay rate T H ^ bb vs. Higgs mass 

In the region of Mh shown in the fig. 3 the three- loop mass corrections are small and does 
not affect the overal picture. However, in the low mass region their impact is significant (fig.4). 

To obtain the decay rate Th^cc one should substitute b — > c everywhere in eq.(|35D except 
the summation bounds, which remain the same. The 6-quark "threshold" crossing effect will 
be represented by the term ~ A(m.&/m c ). 
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fig. 4 The partial decay rate T H 



■ bb 



VS. 



Higgs mass (low mass region) 



Since the derivation of the eqs. (|22|,|34|) has not been restricted by the particular quark flavor, 
one could apply these results to the H — > it decay channel. Q The only restriction that apply is 
that the Higgs mass Mh must be sufficiently larger then 2mt in order the process to be allowed 
and the expansion in terms of mf/Mfj to be a legitimate. Note, that in those equations for 
the it mode one should consider six flavor theory and the summation index v will run over all 
six quark flavors. Thus, from the eq. (0) at mf = m t and Nj = 6 one obtains: 



3y/2G F M H 



H^tt 



-m. 



57T 



Am 2 t 



+ 



a. 



(6) 



{M H ) 



7T 



( a 



(6) 



Hj 



71 



-1.50631 + 247.14204 



m t 



17.41667- 169 



m t 



log 



Ml 

2 



5 (D 



m t 
Ml 



- 2.66667- 32 



m t 
M% 



+ 0.25-27 



mi<m t Km t / m v <m t 1V1 H 



m t 



Mj 



log 



2MI 

2 

mf 



mz 



(37) 



The numerical values of A(mi/m t ) defined in the eq. (|30D are given in the table 1. In the above 
equation the known exact dependence on the top mass at one and two-loop levels is used (see 
eq.(||)). It is straightforward to express the c4 6) in terms of using the relation (0). For 
the conservative lower bound on the top mass m t = 131 GeV and for the Higgs mass at, for 
instance, M H = 400 GeV one obtains: 



The author is grateful to the referee for Physics Letters B for his suggestion. 
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T H ^ ft = — ^-^ m 2 t \ (1 - 0.644 + 0.069 + 0.005 + 0.0008 + ...) 

+ 8 v ; (-1.465 + 4.889 - 1.739 - 0.005 - 0.0003 + ...) 

+ (^£iMA\ 2 (-41.59 + 57.02 + ...)} (38) 

In the above equation at each order of the expansion in rn^/Mjj (at m t = 131 GeV and 
Mjj = 400 GeV) is given. One can see that for the case H —*■ tt the mass corrections are 
important. In fact, the massless limit coefficients (e.g., at 0(a s )) does not reproduce even a 
correct sign. 

6. Scheme-scale dependence and the theoretical uncertainties. The results for 
the decay rates are scheme-scale dependent. The present analysis is restricted by the one 
parametric family of the MS- type schemes |]14j| . The ambiguity is due to the dependence on 



the parameter /1, which enters via dimensional regularization. (For the similar analysis of the 



r 'z^hadrons see [Q.) In the solution of the RG eq.([U|) one replaces fij^g = e t Mj I and a 
where the parametrization of a\ is given by the eq.(E3) with A t = e~ t/2 Ajjo. 



fig. 5 The approximants of the T H ^ b -^ vs the scale parameter t. 

The one-, two- and three- loop approximants for the T H ^ bb in terms of running quark mass 
(eq.(^), with N = 5 and m/ = m&) vs. the scale parameter t are plotted in the fig. 5. One can 
see that the higher order corrections diminish the scale dependence from 40% to nearly 5%. 
The solid curve (corresponding to the three-loop result) became flat in the wide range of the 
logarithmic scale parameter t. Moreover, the choice t — (MS-scheme) satisfies Stevenson's 
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"minimal sensitivity" principle 133]. The mass corrections does not change the overall picture. 
It should be stressed once again that the present consideration is restricted by the MS'-type 
schemes. On the other hand one could carry out the scheme invariant analysis along the lines 
of 11 |35j (see also 0). 



It is known that the dependence on the unphysical scale parameter is due to the truncation 
of perturbation series at particular order. If one sums up to "all orders" , the scale dependence 
will disappear. Thus, one may try to relate the remainder dependence on the scale to the sum 
of all uncalculated terms and estimate the size of the theoretical uncertainty in the reasonably 
wide neighborhood of the initial choice of scale (in our case t = 0). The deviation of the 
approximant from the constant is used as a measure of theoretical error. From the fig. 5 one 
estimates the theoretical uncertainty at 5%. 

The comparison of fig. 3 and fig. 5 shows that the high order corrections resolve the very large 
discrepancy (more then factor 2!) between the results for T H ^ bb in terms of running and pole 
quark masses. The remaining discrepancy is about 5%, which is the same as the estimated 5% 
uncertainty due to the uncalculated higher order terms. 

Besides the above estimated theoretical error, there are two more contributions, with the 
potentially important effects comparable with the present uncertainty. The first most likely 
non negligible contribution comes from the diagrams in fig.2, with the virtual top quark in 
the triangle fermionic loop. (The similar diagrams in the case of Z-boson decay gave sizable 
correction [0.) The other contribution may come from the diagrams in fig.l with the virtual 
top quark. The decoupling theorem may not apply in this case because of possible comparable 
size of the top and the Higgs masses. The contributions from these diagrams (and those in 
fig.2) with virtual bottom and lighter quarks are already included in the result. However, 
in the case of top, one needs to evaluate these diagrams explicitly. The similar diagrams 



for the Z-boson decay gave a somewhat moderate but not entirely negligible correction [|12 
The evaluation of the higher order mass corrections (~ and higher) requires an accurate 
treatment of infrared divergences. On the other hand those corrections are heavily suppressed 
by the powers of m'j/Mf I and numerically will be negligible. The calculation of the 0(a^) 
corrections requires the evaluation of relevant four-loop diagrams up to the finite terms in their 
expansion in Laurent series in e and the method used in to simplify the 0(a^) calculations 
of the decay rate of Z —>■ hadrons does not apply. The above discussion of the theoretical 
uncertainties is mainly for the bb and cc modes. For the case of H — > tt the uncertainty should 
be somewhat higher because of large mass corrections, unless the Higgs mass is in the TeV 
region. 

7. Summary. The quark mass corrections to the decay rate of the SM Higgs particle 
into the quark-antiquark pair to 0(ot 2 s ) perturbative QCD are calculated. The previously 
known 0(a s ) results with explicit mass dependence and the three-loop massless limit result 
with running quark mass parametrization are independently confirmed. The expression for 
the decay rate is obtained both with running and pole quark mass parametrizations. It was 
found that the quark mass corrections are not entirely negligible and they decrease the large 
and negative massless coefficients. For the decay mode H — > tt the mass corrections are large 
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and important. On the other hand, the higher order corrections resolve the large discrepancy 
(more then factor 2) between the results for the decay width in terms of running and pole quark 
masses. It was also found that the three-loop QCD correction reduces the scale dependence 
significantly. The theoretical error of evaluation of the QCD contribution was estimated at 
5% for the decay mode H — > bb plus, the possible effects of the virtual heavy top quark are 
emphasized. 



Finally one should mention the recent work [37|, where the similar problem without calcu- 



lating three-loop mass corrections has been discussed. However, the relation between running 



and pole quark masses used in p7| , seems to be incorrect. 
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